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Here we present a theoretical study of the thermoelectric transport through [2,2]paracyclophane- 
based single-molecule junctions. Combining electronic and vibrational structures, obtained from 
density functional theory (DFT), with nonequilibrium Green’s function techniques, allows us to 
treat both electronic and phononic transport properties at a first-principles level. For the electronic 
part, we include an approximate self-energy correction, based on the DFT+E approach. This 
enables us to make a reliable prediction of all linear response transport coefficients entering the 
thermoelectric figure of merit ZT. Paracyclophane derivatives offer a great flexibility in tuning their 
chemical properties by attaching different functional groups. We show that, for the specific molecule, 
the functional groups mainly influence the thermopower, allowing to tune its sign and absolute 
value. We predict that the functionalization of the bare paracyclophane leads to a largely enhanced 
electronic contribution Z e \T to the figure of merit. Nevertheless, the high phononic contribution to 
the thermal conductance strongly suppresses ZT. Our work demonstrates the importance to include 
the phonon thermal conductance for any realistic estimate of the ZT for off-resonant molecular 
transport junctions. In addition, it shows the possibility of a chemical tuning of the thermoelectric 
properties for a series of available molecules, leading to equally performing hole- and electron¬ 
conducting junctions based on the same molecular framework. 


I. INTRODUCTION 

Nanostructured materials provide a promising 
route towards high-performance thermoelectric energy 
conversion. 1 ^ Nanostructuring does not only allow 
to improve the thermoelectric performance, but also 
enables novel applications like on-chip solutions for 
energy harvesting and refrigeration. 2 ! Recently it became 
feasible to probe the thermoelectric properties and 
energy conversion in molecular junctions.SHE Within the 
different approaches to novel thermoelectric materials, 
molecular junctions are the ultimate level achievable 
regarding device miniaturization and provide a great 
flexibility to tailor their thermoelectric properties, 
e.g. by chemical synthesis. Presently the applicability 
of organic thermoelectric materials is limited by the 
absence of high performance electron-type materials 3 
The possibility to realize equally well performing hole- 
and electron-type (p- and n-type) molecular junctions 
and to enhance their effrciencyEHEIi makes them very 
appealing for further scientific investigation. Irrespective 
of their practical use, they represent ideal systems 
to study thermoelectric transport at the atomic scale 
and to improve the understanding of the fundamental 
processes of nanoscale energy conversion. 

On the theory side the combination of DFT-based elec¬ 
tronic structure calculations with nonequilibrium Green’s 
function techniques allows an atomistic description of 
the electric and electronic thermal transport properties. 
With this approach important trends like the influence 
of the chemical and structural properties on the charge 
transport can be successfully capturedPiUHUl Moreover 
the recently introduced DFT+E method, which improves 


the description of the level alignment at the metal- 
molecule interface, yields a quantitative agreement with 
experimental results at low computational costsP^ 


Recent computational studies on molecular 
junctions 15 '" 1 ' suggest that antiresonances in the 
electronic transmission close to the Fermi energy can 
largely enhance the thermopower. The decrease of 
the electric conductance will be compensated by a 
simultaneous decrease of the electronic contribution to 
the thermal conductance. Consequently, this is reported 
to lead to an increase of the electronic contribution 
to the thermoelectric figure of merit Z e \T. In such 
situations the overall thermoelectric figure of merit ZT 
is eventually limited by the phononic contribution to 
the thermal conductance (see Eqs. (18) and (191 for 
definitions of ZT e \ and ZT). In R.efs. I15UT7| however, 
this phonon thermal conductance is either neglected 
or just taken into account in an approximate manner. 
While it is speculated in Ref. TB] that the phonon 
thermal conductance might just play a minor role, it 
is suggested in Ref. ]T7] that it can have a significant 
contribution to the overall thermal conductance. Based 
on first-principles calculations of the relevant transport 
coefficients, we demonstrate in this work that it is 
indispensable for an accurate description of the energy 
conversion capabilities of molecular junctions to include 
also the phononic contribution. 


As an accurate total energy method, DFT provides the 
ideal tool to study phonon transport from first-principles 
and has already been applied successfully to various bulk 
and nanostructured systemsP^ 20 ! However, for molecular 
junctions first-principles calculations of the phonon ther¬ 
mal conductance have up to now only been performed for 
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interactions.^ As basis set def2-SV(P)^ and the respec¬ 
tive Coulomb fitting basis 35 are used. To ensure an accu¬ 
rate description of the vibrational properties, we use very 
tight convergence criteria. Total energies are converged 
to a precision of better than 10 -9 a.u., while geometry 
optimizations are performed until the change of the max¬ 
imum norm of the Cartesian gradient is below 10 -5 a.u.. 

In the construction of the contact geometries, as shown 
in Fig. []J>, we assume the Au atoms of the electrodes in 
(111) orientation to be fixed at their ideal fee lattice posi¬ 
tions (lattice constant ao = 4.08 A). In contrast, the two 
Au apex atoms and the molecule are fully relaxed. While 
the precise metal-molecule-metal geometry will influence 
the thermoelectric properties, we focus here on the effect 
of the substituents. Due to the time-consuming DFPT 
calculations, we therefore restrict ourselves to one repre¬ 
sentative contact configuration. This contact geometry 
or extended central cluster (ECC)p9 and its division into 
left (L), central (C) and right (R) regions is shown in 
Fig. & for the bare, unsubstituted [2,2]paracyclophane. 


Figure 1. (a) Investigated paracyclophane derivatives with 

X = (H, NH2, OH, NO2, COCF3} in pseudo-para position, (b) 
Example of the used contact geometry, here with X = H. 

a few idealized systems J2H3I Systematic studies for real¬ 
istic single-molecule junctions have yet to be provided. 

We present here a fully first-principles based study 
of the electron and phonon transport properties for a 
series of paracyclophane derivatives connected to gold 
electrodes. We focus on the question, how charge and 
phonon transport can be tuned chemically by intro¬ 
ducing different side-groups, which are either strongly 
electron-withdrawing or electron-donating. The chemical 
structure of the considered [2,2]paracyclophane deriva¬ 
tives is given in Fig. 0b namely [2,2]paracyclophane, 
diamino- [2,2] paracyclophane, dihydroxy- [2,2] paracyclo¬ 
phane, dinitro-[2,2]paracyclophane and ditrifluoroacetyl- 
[2,2]paracyclophane, where we, for simplicity, just ana¬ 
lyze the pseudo-para isomers. For [2,2]paracyclophane 
the single-molecule electrical conductance was measured 
recently.^ 


II. FORMALISM AND METHODS 
A. Electronic, geometric and vibrational structure 


B. Electronic transport 

In the Landauer-Biittiker picture the charge current I 
and the electronic thermal current Q e \ are given by 3 ®2H 

I=2 ^ j ^r el (E) [f(E, Ml , T l ) 

-f(E, /XR, T r )], (1) 

Qd = ^ J dE(E - n)r e i(E)x 

[f(E, Ml, T L ) - f(E, fin, Tr)] . (2) 

Here, t c \(E) is the energy-dependent electron transmis¬ 
sion (e.g. from L through C into R, see Fig. 00 and 
f(E,fj,,T) = {exp[(T — n)/k^T] + 1} _1 is the Fermi 
function, describing the occupation of the left (right) 
electron reservoir at the chemical potential fi l (mr) and 
the temperature Tl (Tr). The energy of the electrons 
is measured relative to the average chemical potential 
fi = (/iL + /ir)/2. This is valid in the linear response 
regime, assuming A/i = /iL — hr and AT = Tl — Tr to 
be infinitesimally small. 

In linear response the thermoelectric transport coeffi¬ 
cients, that is the electrical conductance G , thermopower 
S and electron thermal conductance K e ], can be expressed 
as follow PE o lSl 


The electronic structure and the contact geometries 
are obtained at the DFT level, and the vibrational prop¬ 
erties in the framework of density functional perturba¬ 
tion theory (DFPT). Both procedures are implemented 
in the quantum chemistry software package TURBO- 
MOLE 6 . 4 PHH] We employ the PBE functionaPHU 
with empirical dispersion corrections to the total ground- 
state energy, accounting for long-range van der Waals 


9 e 2 

G = —K 0 , 


S = — 


R-i 

eTK 0 1 


( 3 ) 

( 4 ) 



( 5 ) 
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In the expressions e = |e| is the absolute value of the 
electron charge, h is the Planck constant, k R is the Boltz¬ 
mann constant, and T = (Tp + Tr)/2 is the average 
junction temperature. The coefficients in Eqs. M 
are defined as 

K„ = j d Et a {E) (~mp) (E - fi) n , (6) 

where the chemical potential /i ~ Ep is approximately 
given by the Fermi energy Ep of the Au electrodes. 

The electronic transmission function is expressed 
within the Green’s function formalism) 36 * 37 ! The quasi¬ 
particle energies, entering the electronic Green’s func¬ 
tions, are approximated by the Kohn-Sham eigenvalues, 
including a “self-energy correction”. This self-energy cor¬ 
rection of the DFT+E method^ consists of two parts, 
one that accounts for the underestimation of the elec¬ 
tronic gap of the gas-phase molecule in local and semilo¬ 
cal approximations to DFT, the other one for long-range 
correlations in a junction that stabilize the charge on the 
molecule through image charges in the electrodes. De¬ 
tails regarding the self-energy correction can be found in 
Ref. Ell The electronic contact self-energies, on the other 
hand, are obtained as described in Ref. |3B: To construct 
the electrode surface Green’s function, we use 64 x 64 
fc-points in the transverse direction, which was found to 
be sufficient to obtain converged electronic transport co¬ 
efficients. 


H = 


1 

2 


J 2 pI + 




l 

27? 


£x 


(9) 


Here we have introduced mass-weighted displacement op¬ 
erators q ^ and mass-scaled momentum op¬ 

erators pj = P^t\jMe as conjugate variables, obey¬ 
ing [qt,Px\ = and = \p£,p x \ = o, 

and £ = ( j , c) denotes a Cartesian component c = 
x,y,z of atom j at position Rj = R'- 11 + Qj. The 
phononic system is characterized by its dynamical ma¬ 
trix Ke x = which is the second-order 

mass-weighted derivative (or Hessian) of the total ground 
state energy with respect to the Cartesian atomic coor¬ 
dinates, = 9 2 x E. These harmonic force constants 
are computed within DFPT. 

The local displacement basis allows us to partition the 
dynamical matrix into left (L), central (C), and right (R) 
parts 


( K LL K LC 0 \ 

K = K CL K cc RTcr ■ (10) 

\ 0 K rc -Krr J 

Under the assumption of vanishing coupling between L 
and R, K rl = K lr, = 0, the transm issio n function for 
ballistic phonons can be written a d 43 * 46 * 47 -* 


C. Phononic transport 


r ph (E) = Tr [D* CG (E)A L (E)D a cc (E)A R (E )}, (11) 


Since phonons are chargeless, they can only carry a 
heat current. Similar to the electronic expression, this 
phononic heat current can be calculated in the Landauer- 
Biittiker picturd^®! as 

1 

Q ph = - j d EEr ph (E) [n(E, T L ) - n{E, Tr,)] . (7) 

The corresponding thermal conductance due to the 
phonons is given in linear response by 

= (8 , 

where r p p(E) is the phonon transmission and n(E,T) = 
{exp(E/k R T) — l}” 1 is the Bose function, characterizing 
the phonon reservoirs in the left and right electrodes. 

Similar to the electronic transmission function, those 
of the phonons can also be computed within the 
Green’s function formalism. 37 E3Es] The system of inter¬ 
est here, namely a molecular junction or more gener¬ 
ally a nanoscale conductor, consists of an arbitrary, but 
finite-size scattering region (C), connected to two semi- 
infinite electrodes (L and R). For small displacements 
{Qf} around the atomic equilibrium positions the 

phononic Hamiltonian in the harmonic approximation 
can be written as 


where we have expressed the phonon transmission func¬ 
tion in terms of Green’s functions solely defined on C. 

We use the general definition of the nonequilibrium 
phonon Green’s function on the Keldysh contour 

D Zx(P T ’) = ~i( T cSc {qd T )9x( T ')))’ ( 12 ) 

where the displacement operators q^(r) are in the inter¬ 
action picture, Sc = exp[—i/?i f c drif^r)] is the time 
evolution operator on the Keldysh contour, and Tq is 
the corresponding time-ordering operator) 45 ^ 7 ^ 8 * To cal¬ 
culate the phononic heat current, we need the retarded 
and advanced Green’s functions in energy space. They 
are obtained by a Fourier transform 


1 r°° 

D l’ x ( E ) = yj d tD*£{ty Et / h (13) 

of the real-time functions, derived from Eq. (12) a d 45 * 47 * 


Dbc(E) 


(E + i V ) 2 l-Kcc-ni(E) 


-n t r (e) 


(14) 


with Dbc(E) = Dbc(Ey and an infinitesimal quantity 
r) > 0. For non-interacting phonons in the electrodes 
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the corresponding contact self-energies can be calculated 
exactly as 


n r r (£) = K CY d 1 YY (E)K YC . (15) 

Here, d YY (E) is the surface Green’s function of lead 
Y = {L, R}. Additionally, we have defined the spectral 
density of the electrodes as 


A y = i(ir y - n^) 


(16) 


with n y {e) = ny(B)t. 

Similar to our charge quantum transport approach^!! 
we employ also for the phononic case a cluster-based pro¬ 
cedure. We calculate the dynamical matrix K ecc of a 
large, but finite cluster consisting of the molecule of in¬ 
terest and parts of the electrodes. Using the partitioning 
of the ECC, given in Eq. (10 1 and depicted in Fig. [ljr, 
we extract Kcc, and K r,c from K ecc- What re¬ 

mains to be calculated to determine r p i 1 (i?) are the sur¬ 
face Green’s functions of the electrodes d YY {E). They 
are constructed for perfect semi-infinite crystals. The dy¬ 
namical matrix of the Au electrode is derived from those 
of a spherical cluster, see Fig. UK As in the crystal the 
atoms of the Au cluster are positioned on an fee lattice 
with lattice constant do = 4.08 A. We exploit that for 
large enough clusters the force constants from the cen¬ 
tral atom to its neighbors will be similar to the bulk. 
Following the procedure introduced in Ref. 132] for the 
charge transport, we construct the dynamical matrix of 
the crystalline electrode via these extracted bulk cou¬ 
plings. From this we compute surface and bulk Green’s 
functions, d YY (E) and d Bu]k (E), resp ectively, iteratively 
by means of the decimation technique 1 !| 3S1 

To check the quality of the extracted bulk parameters, 
we compare the bulk phonon density of states (DOS) 
against available experimental data. It can be directly 
obtained from the retarded bulk Green’s function 


p(£) = -Mm(h[<f B , k (E)]). 
7r 


(17) 


The comparison of the DOS, calculated from parame¬ 
ters extracted from Auig, AU 171 and AU 333 clusters, with 
the experimental data of Ref. [50] is shown in Fig. [2 ]d. 
For the smallest cluster Au lg the calculated DOS already 
resembles the general features of the experimental one. 
However, in the lower energy range from 0 to 15 meV 
the agreement is still poor. Increasing the cluster size to 
171 atoms leads to a change especially in the low-energy 
range of the DOS. If we increase the cluster size further 
to 333 atoms, we observe only small modifications. Com¬ 
paring the DOS calculated for the parameters extracted 
from the Au 333 cluster with the experiment,^we see that 
the DFT results resemble the experimental DOS closely, 
but overestimate the vibrational energies by around 15%. 
By scaling the energy axis of the calculated DOS by a 
factor of 0.85, we basically obtain a perfect match with 
the experimental results. Regarding the scaling, we note 




Figure 2. (a) Shape of the three clusters, Auig, Aum and 
AU 333 , used to extract the bulk parameters, (b) Bulk phonon 
density of states calculated with parameters extracted from 
the three cluster shown in (a) and experimental DOS from 
Ref. [Sf7 The blue dashed curve is the DOS as obtained from 
the A 11333 cluster scaled by 0.85 to fit the experimental vibra¬ 
tion energies. 


that in principle both the x- and y-axis of the DOS plot 
need to be scaled simultaneously with inverse factors to 
conserve the energy integral of the DOS. But since the 
experimental data is measured in arbitrary units, a scal¬ 
ing factor for the y axis is already used. 

We want to point out that, in order to obtain a reliable 
DOS as well as phonon transport properties, it is crucial 
to choose the broadening factor 77 sufficiently small. Oth¬ 
erwise the sharp features in the spectra are smeared out. 
Here, we used a broadening of 77 = 6 . 8 /xeV, for which 
we obtained both a converged DOS and converged trans¬ 
mission spectra with respect to 77 . Due to the small y 
a large number of transverse fc-points has to be chosen 
in the decimation procedure. We found the results to be 
converged for a fc-point sampling of 512 x 512. After the 
DFPT step, this large number of fc-points is the bottle 
neck of the transport calculations. However, the calcula¬ 
tion can be parallelized with respect to the ^-points and 
after transforming the surface and bulk Green’s functions 
back to real space, they can be conveniently stored on a 
hard drive and reused in later calculations. 
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Figure 3. Energy-dependent electronic transmission probabil¬ 
ity Tei(E) relative to Ef- 


D. Thermoelectric figure of merit 
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The common measure for thermoelectric efficiency is 
given by the figure of merit 


ZT = 


GS 2 

K e \ + ftph 


T, 


(18) 


which is determined from the thermoelectric transport 
coefficients in Eqs. and ([8]) in the linear response 

regime. An upper bound for ZT in the limit of vanishing 
phonon thermal transport At p h —> 0 is given by the purely 
electronic contribution 


Z el T = —T 

K e \ 


& 

L ' 


(19) 


Here, we have introduced the Lorenz number L = 
K e i/GT. With Z e \T we can express the figure of merit 
in a slightly different form as 


ZT = 


ZelT 

1 "F ^ph/ ^-el 


( 20 ) 


III. RESULTS AND DISCUSSION 

For the unsubstituted paracyclophane molecule 
(X = H), Ef is located almost in the middle of the 
HOMO-LUMO gap (Fig. [3]). The conductance is found 
to be Gdft+s = 0.004 Go (without self-energy correc¬ 
tion Gdft = 0.07Go), which compares reasonably well 
with the experimental data G eX p ~ 0.01 Go, given in 
Ref. [55] Introducing the substituents leads to a shift of 
the molecular orbital energies relative to For the 

electron-donating substituents IMH 2 and OH, the molecu¬ 
lar frontier orbital energies increase due to the increased 
Coulomb repulsion in the conjugated 7r-electron system 
of the benzene rings. This moves the HOMO resonance 


Figure 4. Temperature dependence of the thermoelectric 
transport coefficients, namely (a) electric conductance, (b) 
electron thermal conductance, and (c) thermopower, (d) Elec¬ 
tronic contribution Z e \T to the figure of merit. 


around 1.1 eV closer to Ef- The effect is slightly larger for 
NH 2 than for OH (Fig. [ 3 ]). The LUMO level is largely un¬ 
affected, which suggests that the two substituents mainly 
couple to the occupied 7r-orbitals. For the electron- 
withdrawing substituents NO 2 and COCF 3 , we observe 
the opposite effect. The decreased Coulomb repulsion 
moves the LUMO resonance around 1.5 eV towards Ep, 
with a slightly larger shift for COCF 3 (Fig. |3j. 

In Fig. gjr-c the temperature dependence of the thermo¬ 
electric transport coefficients is summarized. Before dis¬ 
cussing the individual features, it is worth stressing that 
due to the off-resonant transport and smooth transmis¬ 
sion functions around Ep, the temperature dependence 
of the thermoelectric transport coefficients (Eqs. ©-©) 
is dominated by t he low est order of the Sommerfeld ex¬ 
pansion of Eq. (§PE3 

Despite the substituent-induced movement of the or¬ 
bital energies, the transmission at Ep remains roughly 
constant as compared to the unsubstituted molecule. 
Hence the influence of the substituents on the electric 
conductance G as well as the electron thermal conduc¬ 
tance K e \ ss LTG remains small. For NH 2 , G and K e i are 
slightly increased. For the other three substituents a de¬ 
crease is observed (Fig. & ,b). However, the slope of the 
transmission function at Ep , and hence the thermopower 
S, changes substantially. For the unsubstituted molecule 
the transmission function is flat around Ep, resulting in 
an almost vanishing thermopower (Fig. |4jc) . On the other 
hand, for all four substituents the absolute value of S 
is largely increased (Fig. [I]:). Moreover the electron- 


































6 



E [eV] E [eV] E [eV] 


Figure 5. (a) Calculated and measured^] infrared spectra of [2,2]paracyclophane in the gas phase, showing the infrared-active 
molecular vibrations, (b) Phonon DOS of a Au (111) surface, (c) Phonon transmission r p h(-E) for [2,2]paracyclophane connected 
to Au electrodes. Green-shaded areas indicate the energy window, where the electrode phonon DOS is finite. 



Figure 6. Temperature-dependent phonon thermal conduc¬ 
tance. Inset: Ratio Kph/ftei between phonon and electron 
thermal conductance. 


donating substituents NH2 and OH tune the slope of 
r e \(E) at E-p to be negative, giving a positive ther¬ 
mopower, characteristic for HOMO-dominated transport 
(p-type) j 5 l 55 l E3 In contrast, the electron-withdrawing 
substituents give rise to LUMO-dominated transport (n- 
type), and hence the thermopower is negative (Fig. Bo- 

Ultimately, we can combine the thermoelectric trans¬ 
port coefficients to the electronic contribution to the fig¬ 
ure of merit Z e \T , which is displayed in Fig. B- Due 
to the almost vanishing thermopower of the unsubsti¬ 
tuted molecule, Z e \T remains very small. On the other 
hand, due to the quadratic dependence of Z e \T on S, 
the enhanced thermopower for the substituted molecules 
increases Z e \T largely. The N ^-substituted molecule 
shows the highest Z e \T of around 4 • 10 -3 at room tem¬ 
perature (T = 300 K). Overall Z e \T remains rather small 
for the [2,2]paracyclophane derivatives studied here. The 
fact that it is possible to tune the system to be either p- 
or ro-type by simply replacing the substituents is very 
encouraging, as this simple approach can be easily trans¬ 
ferred to other 7r-conducting systems. 


Now we turn to the phononic contribution to the ther¬ 
mal conductance, which is often neglected in the cal¬ 
culations or treated in an approximate manner. Here, 
however, we determine also the phononic thermal con¬ 
ductance accurately from first-principles, as described 
in Sec. |II C Compared to the molecular vibration spec¬ 
trum (Fig. 5^,), the electrode phonon density of states 
(DOS) is finite only in a narrow energy window (Fig. B)- 
Even if the presence of the electrodes can rather strongly 
renormalize the molecular vibration spectrum and the 
character of the individual modes, simply speaking, just 
those molecular vibrations within the energy window de¬ 
fined by the electrode phonon DOS can contribute to the 
phonon transport (Fig. B)- 


The temperature dependence of the phonon thermal 
conductance « p h is displayed in Fig. [6] We obtain the 
expected ballistic behavior with a steep increase for low 
temperatures, as the phonon modes are getting occupied, 
while for higher temperatures K p h approaches a constant 
value. Substituting the paracyclophane has a twofold 
effect. Firstly it influences the molecular vibration spec¬ 
trum directly, moving modes in and out of the energy 
window defined by the electrode phonon DOS. Secondly, 
it changes the binding geometry between the molecule 
and the Au electrode, modifying at the same time the cou¬ 
pling between the molecular vibrations and the electrode 
phonons. Both effects are, however, difficult to separate. 
Substituting with NO 2 increases the phonon thermal con¬ 
ductance, while the other three substituents decrease it 
compared to the unsubstituted molecule. Eventually K p h 
remains comparable for all five molecules (Fig. [6]) and 
is, in the temperature range 10 K < T < 400 K, for all 
molecules at least one order of magnitude larger than the 
electronic thermal conductance (inset of Fig. [6]). 


Since K p h/ft e i > 10 for 10 K < T < 400 K, ZT 
is strongly suppressed according to Eq. (201 (Fig. [7|. 
The largest ZT is obtained for NH 2 substitution with 
ZT = 4 ■ 10 -4 at room temperature (T = 300 K). The 
OH and COCF3 substituted paracyclophanes have al¬ 
most the same ZT in the whole temperature range with 
ZT = 0.9 • 10 -4 at T = 300 K, but the former one is hole¬ 
conducting while the latter one is electron-conducting. 
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Figure 7. Temperature dependence of the figure of merit ZT, 
including both electronic and phononic linear response trans¬ 
port coefficients. 

IV. CONCLUSIONS 

For a series of paracyclophane-based single-molecule 
junctions we computed the electrical and heat trans¬ 
port properties fully from first-principles. Combining 
DFT + S-based electronic structure and DFPT-derived 
harmonic force constants with the Green’s function for¬ 
malism and the Landauer approach allowed us to make 
quantitative predictions of all linear response transport 
coefficients relevant for the thermoelectric figure of merit 
ZT. 

We showed that the transport properties can be tuned 
by functionalizing the paracyclophane molecule appropri¬ 
ately. The main influence of the functional groups was 
found to be on the junction thermopower, where both its 
absolute value and sign are affected. This is due to the 
fact that the Fermi energy lies right in the middle of the 
molecular HOMO-LUMO gap of the bare paracyclophane 
molecule, leading to a vanishing thermopower and con¬ 
sequently ZT. The electron-donating substituents NH 2 
and OH tune the thermopower to be positive, while the 
electron-withdrawing groups I\I0 2 and COCF 3 give rise to 
a negative value. The largest enhancement of the abso¬ 
lute value of the thermopower was found for NH 2 , fol¬ 
lowed by OH and COCF 3 . The latter two substituents 
yield a comparable absolute value of S, but a different 
sign. 

Naturally, the non-vanishing thermopower increases 


the electronic contribution Z e \T to the thermoelectric fig¬ 
ure of merit largely. However, for the studied off-resonant 
conduction, the phonon thermal conductance is, even at 
room temperature, at least one order of magnitude larger 
than the electronic one, leading to a significant suppres¬ 
sion of the overall ZT. Therefore we need to stress that 
it is necessary to include the phonon contribution to the 
thermal conductance for low-conducting molecular junc¬ 
tions in order to obtain accurate and reliable predictions 
of ZT. 

In our present study the largest ZT was found for the 
NH 2 -functionalized paracyclophane with ZT = 4 • 10 -4 
at T = 300 K. Even if the ZT remains rather small, we 
could demonstrate for a series of realistic molecules that 
it is possible to chemically adjust the transport properties 
and to enhance the thermoelectric figure of merit as well 
as to realize p- and n-conducting junctions using the same 
molecular framework. 

In the future more complex molecular structures and 
different functional groups could be used to enhance the 
thermoelectric figure of merit further. Beside the chem¬ 
ical tuning of mainly the electronic transport proper¬ 
ties, which was demonstrated in this work, limiting the 
phonon contribution to the thermal current is necessary 
to improve ZT. Nonlinear effects, leading to increased 
phonon-phonon scattering, could provide a convenient 
way to suppress the phonon transport. Such a discus¬ 
sion will be the topic of future work. 
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